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Preface 


The  purpose  of  this  study  was  to  investigate  a  modeling 
technique  to  simply  and  adequately  model  an  adhesive  bonded 
patch  on  an  orthotropic  plate.  ^  No  previous  models  have  been 
found  which  used  the  following  method:  three-dimensional 
spring  elements  to  model  the  bond. 

This  report  is  limited  to  linear  material  behavior. 

While  this  does  not  technically  tell  the  entire  story  behind 
the  behavior  of  the  patched  plate,  I  believe  that  it  provides 
a  complete  picture  of  the  modeling  process  for,  and 

N 

applicability  to,  actual  composite  repairs^^This  thesis  is 
also  the  basis  for  further  work  into  the  nonlinear  material 
range  and  three-dimensional  stress  states.  ■<‘1^ — 
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Abstract 

-1/ 

To  further  advance  aerospace  battle  damage  repair 
techniques  for  composite  repairs,  a  finite  element  model  of  an 
orthotropic  plate  with  a  hole  and  patch  was  constructed  to 
investigate  the  accuracy  of  such  a  model .  The  fundamental 
problem  in  this  model  is  the  development  of  a  simple  anc 
accurate  method  of  modeling  the  bond  between  patch  and  plate. 
MSC/NASTRAN  spring  elements  were  used  to  this  end.  The  spring 
elements  provided  an  elastic  method  of  load  transfer  between 
the  plate  and  the  patch.  MSC/Nastran  was  used  for  model 
solution  while  SDRC  I-DEA^^was  used  for  Pre-  and  Post¬ 
processing  of  the  finite  element  model.  Two-dimensional 
CQUAD4  elements  were  used  to  model  the  plate  and  patch. 
Physical  test  specimens  were  constructed  to  validate  the 
model .  The  plate  considered  was  a  16  ply  AS4/3501-6 
graph! te/epoxy  composite  with  dimensions  30"  x  10"  x  0.0832" 
with  a  five  inch  diameter  hole  centered  in  the  plate.  Three 
different  patch  concepts  were  modeled.  The  strain  correlation 
between  the  finite  element  model  and  the  test  specimens  was 
excellent  for  points  far  away  from  the  hole  and  centered  on 
the  patch.  Areas  near  the  edge  of  the  patch  and  hole  were  not 
as  successful,  probably  due  to  the  rapidly  changing  strain 
fields  near  these  areas. 
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MODELING  AND  ANALYSIS  OF  ORTHOTROPIC  PLATE 
WITH  CIRCULAR  CUTOUT  AND  ADHESIVE  BONDED  PATCH 

I .  Introduction 

Purpose .  Today's  modern  military  aircraft,  such  as  the  ATF, 
AV~8B,  B-2,  and  ATA,  are  designed  to  very  tight  weight 
tolerances.  The  structure  weight-to-gross  weight  ratio  is 
constantly  decreasing  as  aircraft  design  advances.  These 
advances  are  possible  through  the  use  of  lightweight  laminated 
composite  materials  in  primary  and  secondary  load  bearing 
structures.  The  orthotropic  materials  allow  structure 
material  properties  to  be  oriented  favorably  with  the  load  and 
load  path  the  structure  will  be  required  to  withstand. 
Presently,  this  is  usually  accomplishes  with  graphite/epoxy 
(gr/ep)  structure,  where  the  lay-up  of  the  gr/ep  is  tailored 
to  meet  the  load  requirements. 

These  advances  in  design  and  manufacturing  have  markedly 
increased  the  performance  of  military  aircraft  and  at  the  same 
time  left  the  armed  services  with  a  gap  in  their  battle  damage 
repair  concepts.  The  aircraft  battle  damage  repair  (ABDR) 
concept  stresses  metallic  repairs  for  all  structural  damage 
incurred  on  aircraft  (l;2-3).  This  concept  was  viable  in  the 
era  of  all,  or  mostly,  metallic  aircraft.  Battle  damage  to 
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today's  and  tomorrow's  aircraft  may  require  more  than  just 
aluminum  and  rivets  to  adequately  repair. 

Ob jecti ve .  This  thesis  is  primarily  an  investigation  of 
modeling  techniques  applicable  to  damaged  composite  panels 
with  adhesively  bonded  "scab"  patches.  Noted  below  are 
several  goals  which,  if  accomplished,  could  lead  to  advances 
in  battle  damage  repair  of  composite  structures.  The 
following  goals  of  this  thesis,  in  order  of  importance,  are  as 
foil ows : 

1.  Develop  an  accurate  and  workable  finite  element  model 
(FEM)  of  an  adhesively  bonded  patch  and  plate. 

2.  Characterize  and  evaluate  potential  ABDR  class 
materials  needed  to  accomplish  the  modeled  repairs. 

3.  Evaluate  potential  ABDR  class  hardware  required  to 
accomplish  the  modeled  repairs. 

Procedure .  The  FEM  will  be  "constructed"  using  Structural 
Dynamics  Research  Corporation’s  (SDRC)  Integrated  Design 
Engineering  Analysis  Software  (I-DEAS)  on  a  Digital  Equipment 
Corporation  (DEC)  VAX  11/785  running  the  VMS  operating  system. 
Within  this  environment,  the  FEM  of  the  plate,  patch,  and  bond 
will  be  defined  along  with  the  loads,  restraints,  and 
constraints  of  the  model.  Model  solution  will  be  accomplished 
outside  of  I-DEAS.  I-DEAS  will  be  used  to  write  a  MacNeal- 
Schwendler  Corporation  (MSC)/NASTRAN  data  deck.  The  data  deck 
will  then  be  transferred  to  a  SUN  4  Workstation  running  under 
the  Unix  operating  system.  The  SUN  will  actually  solve  the 
model  using  MSC/NASTRAN.  Once  the  solution  is  complete,  the 
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post-processing  will  take  place  inside  I-DEAS.  The  required 
hard  copy  output  is  then  generated. 

The  physical  model  will  be  manufactured  from  gr/ep  and 
configured  to  match  the  FEM.  Nine  base  panels  will  be 
manufactured.  Eight  will  have  centrally  located  holes 
matching  the  FEM's.  The  panels  are  repaired,  or  left  as  is, 
and  then  strain  gages  are  attached  at  specific  locations.  The 
panels  are  then  tested  in  tension  and  the  strain  gage  data 
recorded . 

The  final  phase  of  the  procedure  involves  reducing  the 
strain  gage  data  and  comparing  it  to  the  data  generated  using 
the  FEM. 
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I I .  Background 


ABDR .  According  to  T  0  l-lH-39,  General  Technical  Manual , 

Aircraft  Battle  Damage  Repair: 

The  primary  purpose  of  ABDR  is  to  restore  sufficient 
strength  and  serviceability  to  permit  damaged  aircraft  to 
fly  additional  operational  sorties,  of  at  least  partial 
mission  capability,  within  time  to  contribute  to  the 
outcome  of  the  on-going  battle.  (l:vii) 

This  manual  is  dedicated  almost  entirely  to  metallic 

structural  repairs,  and  the  portions  dealing  with  damaged 

composite  structures  treat  those  structures  as  basically 

metallic  in  nature,  (i.e.  they  use  the  same  type  of  repairs  as 

for  the  metallic  structures). 

There  are  on-going  efforts  to  determine  the  most 

efficient  method  of  field  level  composite  damage  repair. 

Hinkle  and  Van  Es  (2:11)  recently  completed  a  contract  with 

USAF,  in  which  they  demonstrated  the  accuracy  and  viability  of 

the  McDonnell  Douglas  CREPAIR;  a  computer  based  repair  design 

and  analysis  tool.  This  effort  was  dedicated  entirely  to 

bolted  repairs  of  composite  panels  using  metallic  patches. 

Most  accepted  repair  methods  to  date  still  revolve  around 

metallic  patches  fastened  with  mechanical  bolts.  Adhesively 

bonded  repairs  are  being  investigated  on  the  national  and 

international  level  as  well  (3:6-1  to  6-21). 

There  are  many  factors  to  be  considered  when  preparing  to 

explore  the  patch-bonded-on-orthotropic-plate  concept: 
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1.  Material  analysis  of  the  orthotropic  plate. 

2.  Material  analysis  of  the  adhesive. 

3.  Material  analysis  of  the  patch. 

4.  Surface  preparation  of  patch  and  plate. 

5.  Bonding  environment  requirements. 

6.  Loads  analysis  of  the  plate  with  adhered  patch. 

7.  ABDR  requirements. 

In  the  past  decade  considerable  work  has  been  accomplished  in 
many  of  these  areas. 

Materials .  The  development  of  the  constitutive  matrices  for  a 
composite  laminate  are  well  understood  (4:147-237).  The 
fundamental  building  block  for  an  accurate  composite  material 
analysis  is  the  individual  ply  properties  that  make  up  the 
laminate.  Ply  properties  of  given  materials  are  available 
through  various  sources  (5:335). 

Many  small,  and  even  large,  software  packages  have  been 
developed  to  quickly  analyze  laminates  and  compute  their 
constitutive  matrices  (6).  Using  the  constitutive  matrix, 
the  theory  of  elasticity  can  be  used  to  analyze  relatively 
simple  structural  configurations.  However,  an  FEM  is  needed 
for  structures  with  either  complex  geometry,  boundary 
conditions,  or  loads. 

Information  on  the  material  properties  of  adhesives  in 
various  environments  with  different  forms  of  surface 
preparation  are  available  from  many  sources,  governmental  as 
well  as  private  (7).  Typical  relevant  information  on 
adhesives  documents  the  shear  modulus,  shear  strength,  and 
expected  bond  line  thickness. 
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Modeling ■  The  analysis  of  the  completed  plate  with  bonded 
patch  is  relatively  complex  and  has  been  treated  little  in  the 
past.  In  fact,  the  bond  line  and  the  load  transfer  through 
the  bond  line  is  a  major  modeling  problem.  The  thesis 
completed  by  Smith  (8:35-39)  attempted  to  model  a  gr/ep  patch 
on  a  damaged  aluminum  plate.  His  model  was  fundamentally 
simplified  as  noted  below: 

Inherent  within  both  of  these  repair  models  is 
buried  the  assumption  that  a  perfect  bond  exists  between 
the  aluminum  and  the  repair  patch.  Without  this 
assumption  the  model  would  become  significantly  more 
complex . 

Other  sources  have  also  modeled  the  patch-to-panel  interface 
as  a  perfect  bond  (2:6-1  to  6-21).  In  using  the  perfect  bond 
model,  the  section  of  the  plate  that  is  overlapped  by  the 
patch  is  treated  as  a  single  piece  of  material  and  not  as  two 
separate  entities.  The  material  characteristics  of  both  are 
superimposed  to  form  a  third,  hybrid  material. 

The  modeling  of  the  bond  line  is  of  fundamental 
importance  to  this  current  thesis.  An  FEM  of  the  adhesive  is 
necessary  to  adequately  predict  the  load  transfer  between  the 
patch  and  the  plate  as  well  as  the  failure  load  of  the  bonded 
system.  Once  an  adequate  FEM  has  been  constructed  and 
verified,  various  bonded  repair  concepts  can  be  evaluated 
using  ABDR  criteria. 
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III.  Theory 


Orthotropic  Material  Analysis.  Before  an  FEA  can  proceed,  the 
material  characteristics  of  the  elements  must  be  defined.  An 
orthotropic  material,  which  is  what  the  individual  plys  of  a 
composite  structure  are  assumed  to  be,  has  a  constitutive 
matrix  which  relates  the  six  stresses  to  the  six  strains  in 
the  material  axis  system.  The  six  stresses  are  the  three 
normal  stresses  and  three  shear  stresses;  likewise,  the  six 
strains  are  three  normal  strains  and  three  shear  strains.  An 
orthotropic  constitutive  matrix  contains  nine  independent 
material  constants  which  define  the  behavior  of  the  material . 

o^j  -  Stresses 
-  Strains 

Q^j  -  Reduced  Stiffnesses 

The  reduced  stiffness  matrix  simply  relates  lamina  stresses 
and  strains  and  is  comprised  of  the  engineering  constants, 
i.e.  Young's  modulus,  Poisson's  ratio,  and  shear  modulus.  The 
reduced  stiffness  matrix  is  defined  for  a  particular  lamina 
assuming  the  principal  material  directions  coincide  with  the 
coordinate  system  for  the  structure.  For  a  generally  oriented 
lamina,  the  transformed  reduced  stiffness  matrix  must  be  used. 
This  involves  a  simple  coordinate  transformation  to  the  chosen 
axis  system.  Q-barj^j  designates  the  transformed  matrix.  The 
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next  step  involves  combining  laminae  to  form  a  laminate 
stiffness  matrix. 

For  a  general  laminate  under  plane  stress,  using  the 
Kirchhoff  hypothesis  for  plates,  and  accounting  for  out  of 
plane  displacements  and  curvatures,  the  stress-strain 
relationship  for  the  kth  lamina  in  a  laminate  can  be  defined 
as : 

{o}*=  [5]»{{e»}  +z{K}} 

{e*^}  -  Laminate  Middle  Surface  Strains 
{k}  -  Laminate  Middle  Surface  Curvatures 
The  resultant  forces  and  moments  acting  on  a  laminate  are 
obtained  by  integration  of  the  stresses  in  each  lamina  through 
the  thickness  of  the  laminate. 

{N}  =  {a}  dz  =  f^^iojj^dz 

J-t/2 

{M]  =  {a}  zdz  = 

J-t/2 

{N}  -  Resultant  Laminate  Forces 
{M}  -  Resultant  Laminate  Moments 
t  -  Laminate  Thickness 

-  Distance  From  Midsurface  to  Bottom  of  kth  Lamina 
After  integrating  and  partitioning,  the  stiffness  equations 
may  be  represented  as  follows: 


N 
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where 


a '[ate 


(5) 


Aj^j  -  Extensional  Stiffnesses 
“  Coupling  Stiffnesses 
D^j  -  Bending  Stiffnesses 

The  presence  of  the  Bj^j  term  implies  coupling  between  bending 
and  extension  of  a  laminate.  A  symmetric  laminate,  with 
respect  to  lamina  orientation  about  the  midsurface,  will 
result  in  all  coupling  stiffnesses  being  zero.  This  is  the 
configuration  chosen  for  the  composite  panel  in  this  thesis. 

Finite  Element  Analysis.  The  use  of  finite  element  theory  is 
necessary  in  order  to  analyze  structures  too  complex  to  be 
solved  by  classical  analytical  methods.  The  complexities  can 
arise  from  several  sources: 

1.  Complex  geometry. 

2.  Complex  loads. 
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3.  Complex  boundary  conditions. 

Finite  element  modeling,  as  applied  here,  basically  involves 
the  discretization  of  real  complex  structures  into  a  number  of 
smaller  parts  or  regions  called  elements.  Each  element  is 
connected  to  its  surrounding  elements  at  a  number  of  finite 
points,  called  nodes.  The  only  interaction  of  the  elements  is 
through  the  forces  they  apply  at  the  nodes.  Depending  on  the 
problem  being  modeled  and  the  elements  used,  the  nodes  may 
have  from  one  to  six  degrees  of  freedom  (DOF).  The  assumption 
of  small  displacements  and  small  rotations,  and  linear 
material  behavior  lead  to  a  linearized  formulation  and  a 
solution  through  linear  algebra. 

Cook  breaks  down  the  finite  element  analysis  (FEA)  into 
seven  parts: 

1.  Divide  the  structure  into  finite  elements. 

2.  Formulate  the  properties  of  each  element.  This  means 
determining  nodal  loads  associated  with  all  element 
deformation  states  that  are  allowed. 

3.  Assemble  elements  to  obtain  the  finite  element  of  the 
structure . 

4.  Apply  the  known  loads:  nodal  forces  and/or  moments. 

5.  Specify  how  the  structure  is  supported.  This  step 
involves  setting  several  nodal  displacements  to  known 
values . 

6.  Solve  the  simultaneous  linear  algebraic  equations  to 
determine  nodal  DOF  (nodal  displacements  in  stress 
analysis) 

7.  In  stress  analysis,  calculate  element  strains  from 
the  nodal  DOF  and  the  element  displacement  field 
interpolation,  and  finally  calculate  stresses  from 
strains.  (9:4) 
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In  a  FEA,  solutions  are  calculated  for  the  nodal 
displacements  and  rotations  in  all  the  elements  making  up  the 
structure,  for  a  given  loading.  The  elements  making  up  the 
model  must  satisfy  the  conditions  of  equilibrium  and 
compatibility.  That  is  to  say,  the  summation  of  forces  at  any 
node  must  be  zero  and  the  displacements  of  the  nodes  must  be 
such  that  voids  or  displacement  discontinuities  do  not  occur. 

As  referred  to  above,  each  type  of  element  has  its  own 
elemental  stiffness  matrix.  This  matrix  relates  displacements 
for  the  given  DOF  of  the  element  to  the  force  required  to 
achieve  that  displacement.  Once  an  element  has  been  chosen 
that  best  models  the  physical  situation  at  hand,  the 
"stiffness  matrix"  for  that  element  can  be  derived  by  a  number 
of  methods.  Some  common  methods  are  direct  stiffness,  virtual 
work,  and  isoparametric.  The  matrix  equation  relating 
structure  stiffness  ,  nodal  displacement,  and  nodal  loads  is 
shown  below: 

[k:\{d)  =  {r)  (9) 

[K]  -  Structure  Stiffness  Matrix 

{D}  -  Nodal  DOF  Displacement  Vector 

{R}  -  Applied  and  Reaction  Nodal  Loads 

The  physical  meaning  of  [K],  as  well  as  a  procedure  for 

formulating  [K],  is  described  in  the  following  statement: 

The  jth  column  of  [K]  is  the  vector  of  loads  that  must  be 
applied  to  the  nodal  DOF  in  order  to  maintain  the 
deformation  state  associated  with  the  unit  value  of  DOF  j 
while  all  other  nodal  DOF  are  zero.  (9:34) 


The  above  description  for  evaluating  [K]  is  the  direct 
stiffness  approach.  This  approach  is  useful  for  developing  an 
intuitive  understanding  of  [K]  and  evaluating  simple 
structures.  Note  that  [K]  is  singular  in  its  original  form. 

In  order  to  solve  for  unknown  reactions  and  displacements 
Equation  9  must  be  partitioned.  After  known  nodal  DOF 
displacements  and  loads  are  identified,  Equation  9  is 
partitioned  according  to  the  known  and  unknown  portions  of  the 
nodal  displacements  and  forces: 


^11 

^12 

feel 

^2. 

Pcj 

1 

{D^}  -  Known  Nodal  Displacements 
{R^}  -  Applied  Nodal  Forces 
{Djj}  -  Unknown  Nodal  Displacements 
{Rjj}  -  Nodal  Reaction  Forces 

The  above  equation  can  be  successfully  solved  provided  that 
there  are  sufficient  boundary  conditions  to  prevent  rigid  body 
motion  of  the  modeled  structure;  i.e.  determines  if  a 

solution  is  possible.  Now,  solutions  are  calculated  for  the 
unknown  displacements: 

{DJ  =  [Jqj  (  {R^}  -  {DJ  )  (11) 

And  similarly,  the  nodal  reaction  forces  are  calculated: 
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{J?,}  =  [Ji^j  {D,)  *  {Z3„) 


(12) 


Generalized  Stiffness  Matrix.  Deriving  the  elemental 
stiffness  matrix  for  a  given  element  type  is  the  fundamental 
theoretical  procedure  behind  a  PEA.  Actual  computer  program 
equations  start  with  the  elemental  stiffness  matrix  as  a 
known,  dependent  only  on  the  type  of  element  modeled.  The 
potential  energy  form  of  solution  is  a  classical  method. 

lCp=£7+Q  (13) 

Xp  -  Potential  Energy  on  an  Elastic  Body 
U  -  Strain  Energy 
Q  -  Potential  of  Loads 

For  linearly  elastic  body  under  conservative  loads,  with 
volume  V  and  surface  area  S: 

!!_=/■  {e)dy-f  {u}^{f’}dV 

•’V'2  Jv  (14) 

-/^{u}^{4>}dS-{D}’'{P} 

[E]  -  Constitutive  Matrix 
{6}  -  Strain  Field  Vector 
{F}  -  Vector  of  Body  Forces 
{§}  -  Vector  of  Surface  Tractions 
{P}  -  Vector  of  Loads  Applied  to  DOF 
(u)  -  Displacement  Field  Vector 
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NOTE:  The  above  equation  neglects  initial  strains  and 

stresses . 

The  principle  of  stationary  potential  energy  states  that 
equilibrium  prevails  when  65tp=0  for  any  small  admissible 
variation  of  the  original  structure.  Jip  is  a  function  of  the 
displacement  DOF  defined  by  the  structure.  Applying  the 
principle  of  stationary  potential  energy: 


•  ^  n 


i=l,2,  ...  ,12 


(15) 


-  ith  Structure  Displacement  DOF 
If  Dj^  is  a  nodal  displacement  (or  rotation),  the  above 
equation  is  a  nodal  equilibrium  equation  stating  that  forces 
(or  moments)  applied  at  the  node  sum  to  zero  in  the  direction 
Dj^.  The  potential  energy  equation  components,  for  a  single 
element,  can  also  be  expressed  as: 


U=^{D}'^[FC]  {D}  ;  Q  =  -{D}^{R}  (16) 


This  leads  to  the  equation  for  [K],  in  terms  of  Dj^: 


dDfiDj 


(17) 


This  provides  the  basis  for  generating  the  stiffness 
matrix  for  a  generalized  element.  Typically  though,  the  finite 
element  form  of  the  Rayl eigh-Ritz  method  uses  generalized 
coordinates  which  provide  for  a  simpler  mathematical 
evaluation  of  the  stiffness  matrix: 
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1  m  n 

2=1  2=^+1  i=in+l 

u,  V,  and  w  -  Displacement  Components  of  a  Point 
aj^  -  Generalized  Coordinates 
=  fj^{x,y,z) 

Each  of  the  functions  fj^  must  be  admissible,  i.e.  each  must 
satisfy  compatibility  and  essential  boundary  conditions. 
Usually  the  fj^  are  polynomials.  Now  becomes  a  function  of 
DOF  and  the  equilibrium  equation  is  defined  by  n  algebraic 
equations : 

du- 

-=-£=0  for  2=1, 2, 3  /  .  .  . /^  (19) 

3a  i 

A  Ray  1 eigh-Ritz  solution  is  usually  approximate  because  the 
functions  f^^  are  usually  incapable  of  exactly  representing  the 
actual  displacements.  Additional  sources  of  information  on 
FEA  are  available  (9,  10). 
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IV .  Model inq 


In  order  to  adequately  design  and  test  the  bonded  patch 
repair  concept,  an  analytical  model  must  be  developed.  An  FEM 
which  can  model  the  plate,  patch,  and  the  bond  joining  them 
both  is  necessary.  The  model  must  be  relatively  simple 
(capable  of  being  formulated  and  solved  in  a  reasonable  amount 
of  time  with  computers  at  hand)  and  still  achieve  reasonable 
accuracy  compared  to  the  physical  situation.  The  results  of 
the  model  solution  will  be  verified  through  testing  of 
physical  specimens. 

Physical  Model .  The  plate  configuration  (Figure  1)  was  chosen 
for  ease  of  manufacturing  and  testing,  for  both  the  analytical 
and  physical  models,  as  well  as  being  relevant  to  the 
discussion  of  ABDR  of  composite  structures.  The  dimensions  of 
the  composite  panel  reflect  the  current  limitations  of 
approved  ABDR  composite  repair  procedures.  The  circular 
cutout  is  five  inches  in  diameter  and  centered  on  the  panel . 

A  three-to-one  ratio  of  1 ength- to-width  of  the  panel  is 
considered  fairly  standard  in  order  to  limit  end  effects  from 
the  gripped  end.  The  extra  10  inches  in  panel  length  is  taken 
up  in  the  grips. 

The  material  composition  of  the  panel  is  representative 
of  current  composite  structures;  Hercules  AS4/3501-6  is  a  very 
common  prepreg  used  in  manufacturing  composite  components. 
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All  Dicnensions  are 
in  Inches 


Figure  1.  Composite  Panel  Configuration 
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The  ply  lay-up  of  the  panels  is  (±45,  O2,  90,  O2/  90)^  (Figure 
2),  simulating  a  strength  critical  composite  panel.  Nine 
panels  were  cured  according  to  the  manufacturers 
specifications  (5). 

The  test  matrix  for  the  panels  is  noted  in  Table  1.  The 
undamaged  panel  and  the  two  panels  without  repairs  are  used  to 
verify  the  analytical  model  before  the  patch  and  bondline 
elements  are  added.  If  the  FEM  does  not  accurately  predict 
the  strain  fields  of  these  relatively  simple  cases,  then  the 
FEM  with  the  patch  and  the  bondline  included  could  not  be 
expected  to  accurately  model  the  physical  specimen. 


Number  of  Panels 

Repair  Concept 

1 

Undamaged 

2 

No  Repair 

2 

4-Ply  Patch,  7"  Dia 

2 

6-Ply  Patch,  7"  Dia  | 

1  2 

8-Ply  Patch,  9"  Dia  | 

Table  1.  Number  and  Type  of  Test  Panels 


Originally,  the  patch  material  was  to  be  fabricated  from 
3M's  SP  377  CUl,  a  low  energy  cure  unidirectional  gr/ep 
prepreg  tape  (7).  This  material  is  a  possible  ABDR  candidate 
composite  repair  material  due  to  its  relatively  low 
temperature  requirements  for  cure.  Not  only  can  the  material 
be  fully  cured  at  205  ®F  for  two  hours,  but  it  has  a  room 
temperature  storable  life  of  one  year.  The  adhesive  used  to 
bond  the  patch  and  panel  was  also  from  this  family  of  low 
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AS4/3S01-6  TEST  LAM  NATE 

« 

■  0  4 


A  •  Matrix  (Uanbrina) 


1.033t4Ef06 

1.22804Et05 

0.00105590 

1.223D4Et05 

6.27324Et05 

0.0179341 

0.00105590 

0.0179341 

1.73371E+05 

B  -  Mitri 

X  (Mantrane  • 

Bandi ng) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

D  • 

Matrix  (Bending) 

538. 672 

139.110 

19.2031 

139.110 

282.  630 

1^.2091 

19.2031 

19.  2031 

168.280 

S  •  Matrix  (Transvaria  Shaar) 

0.0  0.0 

0.0  0.0 


Larri  nata  thicknaas  -  0.0332  ii 


Figure  2.  Composite  Panel  Properties 
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energy  cure  prepregs.  The  adhesive,  AF  377,  is  an  epoxy  film 
which  is  also  cured  at  205  «F  for  two  hours.  The  film 
incorporates  an  extremely  light  weight  scrim  cloth  to  hold  the 
uncured  epoxy  together.  The  scrim  cloth  also  has  a  secondary 
function  of  controlling  the  bondline  thickness  of  the  cured 
adhesive . 

The  patch  laminates  are  (0)^,  where  X  indicates  the 
number  of  lamina  in  the  patch  (Figures  3,4,  and  5).  The 
patches  are  intended  to  bring  back  the  stiffness  in  the 
longitudinal  direction  only.  Also,  the  patches  are  produced 
such  that  all  the  lamina  of  the  patch  are  the  same  size,  i.e. 
no  taper  of  the  patch.  The  patches  are  centered  on  the 
"damaged"  area. 

The  methodology  of  the  "repair"  is  as  follows: 

1.  The  panel  has  a  five  inch  circular  cutout  centered  on 
it.  The  cutout  was  made  first  with  a  hole-saw  and  taken 
to  the  final  dimensions  with  a  router. 

2  The  patch  plys  are  cut  from  the  patch  material,  SP  377 
CUl . 

3.  The  plys  are  stacked  as  required  and  cured  under 
vacuum,  27  inches  of  Mercury,  and  heat.  The  process  is 
controlled  by  a  prototype  hot-bonder.  This  foot-locker 
sized  machine  controls  the  vacuum  and  heat  applied  to  the 
patch  to  match  a  pre-programmed  cycle  (Figure  6). 

4.  The  panel  is  cleaned  and  abraded  around  the  cutout  in 
preparation  for  bonding  the  patch. 

5.  The  adhesive  film  is  cut  from  a  roll  of  AF  377  and 
temporarily  spot  bonded  to  the  patch. 

6.  The  patch  and  adhesive  film  is  positioned  on  the 
panel  and  taped  in  place.  A  filler  material  is 
positioned  in  the  cutout  to  prevent  the  patch  from  being 
pushed  down  into  the  hole  once  the  vacuum  is  applied. 
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PATCH  4-  PLY 


A  •  Mil  fix  (Ibntrane) 


4.33007Ef05  8158.11 

0.0 

8158.11 

27193.7 

0.0 

0.0 

0.0 

20800.0 

B  • 

Mitrix  (Mentrane  • 

Bendi ng) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

D  ■  Mitrix  (Bending) 

15.61U  0.  294127  0.0 

0.294127  0.980424  0.0 

0.0  0.0  0.  749909 

S  -  Mitrix  (Transvirie  Shear] 
17339.3  0.0 

0.0  17333.3 


Laniftite  thickness  -  0.0208  in 


Figure  3.  Patch  Properties  4-Ply 
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PATCH  6- PLY 


I 

•  0  + 


A  - 

Matrix  (Mambrane) 

8.49511E405 

12237.2 

0.0 

12237.2 

40790.6 

0.0 

0.0 

0.0 

31200.0 

B  •  Matrix  (Mantrane  • 

Bondi  ng) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

D  • 

Matrix  (Bonding) 

52.6884 

0.992679 

0.0 

0.992679 

3.30893 

0.0 

0.0 

0.0 

2.53094 

S  •  MBtrix  (Trinsverse  Shear) 
28000.0  0.0 

0.0  26000.0 


i  n 


Figure  4.  Patch  Properties  6-Ply 
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PATCH  8-  PLY 


I 

.  0  ♦ 


LiminaU  thickness  ■  0.0416  in 


A  -  Mitri X  (Msibrane) 


e.6601SEt05 

16316.2 

0.0 

16316.2 

54387.4 

0.0 

0.0 

0.0 

41600.0 

B  • 

Matrix  (Mantrane  • 

Bendi  ng) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

D  • 

Matrix  (Bending) 

124.691 

2.35302 

0.0 

2.35302 

7.64339 

0.0 

0.0 

0.0 

5. 99928 

S  •  Matrix  (Transvaria  Shear) 
34666.7  0.0 

0.0  34666.7 


Figure  5.  Patch  Properties,  8-Ply 
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7.  The  secondary  adhesive  bonding  is  accomplished  in  the 

same  manner  as  the  cure  of  the  patch. 

Once  the  patch  is  bonded  to  the  panel,  the  conglomeration  is 
ready  for  instrumentation  and  testing. 

The  above  sequence  of  events  was  the  original  plan,  using 
actual  ABDR  capable  material  and  equipment  in  fabricating  the 
composite  "repair."  Unfortunately  this  plan  went  awry.  After 
several  patches  were  cured,  they  were  inspected  by  ultrasonic 
C-scan,  using  a  SRL  Model  1750  Mini-C,  in  a  water  immersion 
test.  The  C-scan  showed  regular  areas  (dark  areas)  of  high 
void  content  in  all  of  the  patches  (Figure  7).  An  inspection 
of  cross-sections  of  one  of  the  patches  revealed  that  the 
prepreg  material  itself  was  resin  starved.  The  cross- 
sectioned  patch  delaminated  within  one  of  the  plys  and  not  at 
a  ply  interface  (Figure  8).  The  SP  377  CUl  prepreg,  which  had 
been  delivered  from  the  manufacturer  about  three  weeks  prior 
to  this  finding,  had  visible  (once  we  knew  what  to  look  for) 
areas  of  resin  starvation  all  along  the  roll  of  material. 

This  problem  made  the  material  useless  for  this  project. 

Due  to  time  constraints,  the  patch  material  used  in  all 
the  patches  was  of  the  same  material  as  the  original  panels, 
AS4/3501-6.  The  configuration  of  the  patches  did  not  change. 
Although  the  material  used  in  adhering  the  patch  to  the  panel 
was  still  AF  377,  the  method  of  bonding  changed.  Instead  of 
the  hot-bonder,  a  hot  press  was  used  for  the  bonding  process 
(Figure  9).  This  was  due  to  time  constraints  and  probability 
of  success. 
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Figure  8.  Micrograph  of  Cross-Section  of  4-Ply  Patch 


The  material  properties  of  AS4/3501-6,  as  obtained 
through  Mall  (11),  are  listed  in  Table  2.  His  data  is  used 
rather  than  the  manufacturers  data  on  the  basis  of 
conservatism  and  completeness. 
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Analytical  Model.  The  analytical  model  is  constructed  with 
the  use  of  SDRC  I-DEAS.  This  software  package  enables  pre- 
and  post -processing  of  the  analytical  model  in  a  graphical 
environment.  I-DEAS  is  used  to  generate  the  Bulk  Data  Deck 
and  Case  Control  Deck  of  a  MSC/NASTRAN  FEM.  The  user  supplies 
the  Executive  Control  Deck.  The  solution  is  returned  to  I- 
DEAS  for  post-processing  and  evaluation. 

MSC/NASTRAN  Version  65  was  used  instead  of  Version  66  as 
the  latter  version,  apparently,  cannot  solve  for  the  inter- 
la^'inar  stresses  or  strains  of  a  laminate  material.  A  data 


- - 

Modulus  of  Elasticity 

El 

20.7E06  psi 

Modulus  of  Elasticity 

E2  and  E3 

1.3E06  psi 

Possion's  Ratio 

V12,  V23*,  and  Vi3* 

0.3 

Shear  Modulus 

^12'  *^23  '  *^13'*^ 

. . . . . 

1.0E06  psi 

Table  2.  Ply  Properties  of  AS4/3601-6 
*  -  Assumed  Property 
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deck  submitted  under  Version  65,  using  Alter  79  for  Rigid 
Format  24,  produces  the  required  results;  while,  the  same  data 
deck  submitted  under  Version  66  will  not  successfully  finish. 

The  first  step  in  the  analytical  modeling  process  is  the 
definition  of  the  panel.  Due  to  the  loading  and  physical 
symmetry,  only  a  quarter  of  the  panel  must  be  modeled.  This 
saves  modeling  and  computation  time.  The  correct  geometry  of 
the  quarter  panel  is  described  as  a  surface.  This  surface  is 
subdivided  into  the  many  finite  elements  that  will  m^del  the 
panel  (Figure  10).  If  bending  moments  had  been  applied,  the 
argument  of  symmetry  would  not  be  valid  due  to  the  ±45  plys  in 
the  panel . 

Once  the  base  panel  is  meshed,  the  element  type  to  be 
used  must  be  defined.  The  element  used  in  the  model  of  the 
panel  is  the  MSC/NASTRAN  CQUAD4  element.  The  CQUAD4  is  an 
isoparametric  membrane-bending  element.  The  forces,  stresses, 
and  the  associated  coordinate  system  supported  by  this  element 
are  illustrated  in  Figure  11  (12:7-9).  In  this  model  the 
membrane  and  bending  properties  are  not  coupled.  This  is 
basically  a  two-dimensional  element.  This  model  is  designed 
to  keep  all  relevant  DOF  in-plane,  i.e.  bending  is  not 
present . 

The  material  properties  associated  with  the  CQUAD4 
elements  are  defined  using  the  material  properties  of  one  ply 
of  AS4/3501-6.  I-DEAS  uses  the  data  on  one  ply  of  material  as 
the  building  block  for  the  material  property  of  an  entire 
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Figure  10.  FEM  Mesh  of  the  Composite  Panel 
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laminate.  The  laminate  is  defined  by  the  number  of  plys, 
orientation  of  the  plys,  and  stacking  sequence. 

After  the  elements  forming  the  panel  have  been  described, 
the  restraints  on  the  model  must  be  applied.  The  restraints 
are  imposed  along  the  two  lines  of  symmetry.  They  keep  the 
model  from  displacing  in  a  rigid-body  fashion  and  force  the 
symmetry  that  would  exist  in  a  model  of  the  whole  panel 
(Figure  12 ) . 

The  loads  are  applied  in  a  consistent  nodal  approach  that 
accurately  reproduces  the  distributed  load  of  a  grip  (on  the 
physical  specimen)  onto  the  boundary  nodes  of  the  model.  The 
loads  are  shown  in  Figure  13.  The  nodal  loads  modeled 
represent  an  equivalent  end  loading,  for  an  entire  panel,  of 
1800  pounds. 

Once  the  restraints  and  loads  have  been  applied, 
constraints  must  be  defined.  If  no  constraints  are  applied, 
the  deformed  model  will  have  an  elongated  central  portion  at 
the  end  of  the  panel  due  to  the  nodal  loads.  This  obviously 
is  not  what  happens  in  the  actual  panel  and  the  constraints 
simply  force  the  end  of  the  model  to  behave  in  the  same  manner 
as  the  physical  panel .  Thus  the  end  of  the  panel  remains 
square  after  loading  (Figure  14).  The  encircled  "MP" 
identifies  the  illustrated  constraint  as  an  MSC/NASTRAN 
multipoint  constraint. 

I -DBAS  is  then  used  to  produce  the  MSC/NASTRAN  Data  Deck 
and  the  solution  is  accomplished  on  a  SUN  4  Workstation.  The 
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Figure  13.  Nodal  Loads  of  the  Panel 
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Figure  14.  Quarter-Panel  Constraints 
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workstation  was  used  for  the  MSC/NASTRAN  solution  instead  of 
the  VAX  as  the  VAX  is  fairly  slow  to  solve  the  computational 
problems  due  to  the  number  of  users,  and  various  demands  on 
the  system.  Afterwards,  data-loader  software  is  used  to 
reformat  the  MSC/NASTRAN  solution  into  the  I-DEAS  universal 
file  format.  The  solution  universal  file  is  then  read  into  I- 
DEAS  for  the  post-processing  phase. 

Before  proceeding  further,  the  model  must  be  tested  for 
convergence  and  accuracy.  The  theoretical  value  for  the 
maximum  stress  of  a  plate  with  hole  under  axial  loading  is 
available  for  homogenous  materials  (13:81).  The  material 
properties  for  the  plate  are  temporarily  redefined  for  those 
of  aluminum.  The  solution  of  the  maximum  stress  on  the 
homogenous  plate  is  compared  to  the  theoretical  value.  The 
FEM  solution  obtains  a  value  of  maximum  stress  within  six 
percent  of  the  theoretical  value.  Next,  the  model  is  tested 
for  convergence;  that  is,  with  a  finer  grid  and  more  elements 
will  the  solution  converge  toward  the  correct  answer.  Figure 
15  is  the  refined  model .  The  maximum  stress  value  for  this 
model  was  within  three  percent  of  the  theoretical  value. 

The  FEM  of  the  patch  is  accomplished  in  the  same  manner 
as  the  panel.  The  patch  FEM  approximates  a  quarter  circle  as 
much  as  the  elements  will  allow;  Figure  16  represents  the 
seven  inch  diameter  patch.  The  patch  is  offset  above  the 
plane  of  the  panel  by  a  distance  equivalent  to  the  adhesive 
thickness,  0.004  inches.  This  is  not  actually  necessary,  but 
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it  eases  the  modeling  and  interpretation  of  the  results  in  a 
graphical  environment.  The  elements  and  nodes  of  the  patch 
FEM  are  modeled  to  match  the  panel  FEM  elements  and  nodes. 
Nodes  of  the  patch  lie  directly  above  nodes  of  the  panel 
(Figure  17).  The  model  is  now  ready  for  the  adhesive 
elements . 

Modeling  the  Adhesive.  The  adhesive  is  modeled  with  the 
MSC/NASTRAN  CELAS2  scalar  spring  element.  This  element  was 
chosen  as  it  can  be  used  as  the  basic  building  block  for  an 
elastic  material  characterization  FEM.  The  basic  nature  of 
the  spring  element,  and  its  parameters,  allows  fine  control 
over  the  FEM  incorporating  it. 

The  element  functions  as  a  linear  spring  associated  with 
two  discrete  points  (nodes).  Any  CELAS2  element  can  act  along 
any  one  of  the  six  DOF  available  to  the  nodes.  That  is,  a 
CELAS2  element  can  be  defined  (between  two  nodes)  to  function 
only  in  the  x-axis  direction;  any  relative  movement  between 
the  two  nodes  in  the  y-  and  z-axis  directions  will  not 
generate  any  reactions  due  to  that  movement.  Therefore,  three 
displacement,  CELAS2  elements  can  be  associated  with  a  pair  of 
nodes.  Each  of  the  three  elements  acts  along  a  different 
axis,  and  the  stiffness  of  any  element  is  independent  of  the 
other  two  (Figure  18).  The  "spring"  between  two  nodes  can 
thus  be  tailored  to  fit  the  behavior  desired. 

The  underlying  assumption  in  using  the  CELAS2  elements  is 
that  the  adhesive  stiffness  is  constant.  This  limits  the  use 
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Kx,  Ky,  AND  Kz 
ARE  INDEPENDENT 
VALUES  COMPRISING 
K 


Figure  18.  CELAS2  Elements  Associated  With  Two  Nodes 
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of  this  method  to  relatively  linearly  elastic  adhesives  below 

the  yield  point.  To  complete  the  model,  scalar  values  of 

stiffness  associated  with  the  spring  elements  must  be 

determined.  Both  these  concerns  will  be  addressed  with  a 

standardized  adhesive  test. 

Calibration  Test.  Testing  of  bulk  adhesives  is 

relatively  straightforward,  with  results  for  many  adhesives 

available  (5).  At  the  time  of  this  thesis  no  neat  resin  data 

was  available  for  3M's  AF  377.  If  it  had  been,  it  would  have 

been  of  marginal  use.  According  to  Hughes,  et  al . : 

Neat  specimens  are  fabricated  in  the  typical  "dog-bone" 
shape  used  for  metallic  materials.  They  are  tested  in 
uniaxial  tension  at  several  strain  rates  and 
temperatures.  This  information,  however,  is  not  directly 
applicable  to  bonded  joint  analysis.  The  constraint  that 
the  adherend  applies  to  the  adhesive  as  it  is  strained 
prevents  the  normal  Poisson  effect  from  taking  place. 
(14:141-175) 

Since  no  adhesive  shear  data  is  available,  at  the  time  of  this 
thesis,  the  only  way  to  gather  the  necessary  information  for 
assigning  values  to  the  spring  elements  is  through  actual 
shear  tests. 

In  order  to  "calibrate"  the  springs,  an  FEM  of  a 
ANSI/ASTM  D  3528-76  test  is  created  using  the  same  methodology 
as  the  patch  and  panel  before.  This  test  is  the  "Standard 
Test  Method  for  Strength  Properties  of  Double  Lap  Shear 
Adhesive  Joints  by  Tension  Loading"  (15:331-336).  Figure  19 
represents  the  specimen  used  in  this  test.  The  FEM  analog  is 
shown  in  Figures  20,  21,  and  22.  The  FEM  of  the  double  lap 
shear  specimen  is  such  that  all  the  elements  have  the  same 
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area  (this  simplifies  the  calibration  effort).  The  nodes  of 
the  FEM  are  defined  such  that  linked  pairs  have  the  same  x  and 
y  coordinate  but  are  separated  in  the  z  (thickness)  direction 
by  the  thickness  of  the  bond,  as  is  the  patch  and  panel.  Each 
pair  of  corresponding  nodes,  of  the  thick  and  thin  adherends, 
are  "connected"  by  three  linear  spring  elements,  as  described 
previously,  though  only  one  of  the  three  is  represented  in 


Figure  22. 


Considering  the  adhesive  thickness,  the  equivalent  spring 
out-of-plane  stiffness  for  a  film  adhesive  is  on  the  order  of 
60%  of  the  equivalent  spring  in-plane  stiffness  (16:160). 

Both  in-plane  springs  have  the  same  value  of  stiffness,  since 
the  adhesive  behaves  similarly  in  both  in-plane  directions. 

The  FEM  spring  values  are  assigned  according  to  the  bond  area 
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Figure  19.  Form  and  Dimensions  of  Double  Lap  Shear  Specimen 

associated  with  the  node  in  question.  Thus,  interior  elements 
have  a  value  of  K  for  the  spring  stiffness,  corner  elements 
have  1/4  K  and  other  boundary  elements  have  1/2  K. 

The  spring  constant  is  defined  using  test  data  of  the 
relative  displacement  of  two  points  on  the  specimen.  The 
points  are  at  either  end  of  the  double  lap,  one  on  the  thick 


Figure  21.  Magnification  of  FEM  Double  Lap  Shear  Specimen 

adherend  and  the  other  on  a  thin  adherend.  The  data  received 
from  the  test  are  strain  and  load.  The  strain  is  converted  to 
a  displacement,  utilizing  the  known  reference  length  of  the 
extensometer  (which  is  the  joint  length).  The  model  is  loaded 
to  the  test  load  (using  a  consistent  nodal  approach  and  proper 
constraints  and  restraints)  and  the  spring  values  of  the  model 
are  adjusted  until  the  displacement  of  the  two  points  in  the 
model  matches  those  of  the  test  specimen.  Once  the  model  and 
specimen  agree,  the  in-plane  spring  value  for  interior  nodes 
is  calculated  for  one  square  inch  of  adhered  area. 

The  spring  stiffness  value  for  one  square  inch  of 
adhesive  is  then  used  to  calculate  stiffness  values  for  all 
the  modeled  springs  based  on  the  bond  area  associated  with  a 
given  pair  of  linked  nodes.  The  bond  area  associated  with  the 
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Figure  22. 


Spring  Elements  of  Double  Lap  PEM 


boundary  nodes  of  the  patch  is  the  equivalent  area  of  the 
actual  patch  and  not  the  elements  describing  the  patch.  Thus, 
although  the  FEM  patch  is  not  round  like  the  physical  patch, 
the  modeled  adhesive  is  representative  of  the  physical 
adhesive  area. 
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V .  Testing 


The  testing  involved  in  this  thesis  is  necessary  to 
determine  the  characteristics  of  the  epoxy  film,  which 
includes  the  linear  or  non-linear  behavior  of  the  adhesive. 
Also  the  testing  will  verify  the  panel  and  patch  FEM  using 
MSC/NASTRAN  CELAS2  spring  elements  to  model  the  bond  between 
them. 

Double  Lap  Shear  Test.  The  CELAS2  elements  must  be 
"calibrated"  before  the  modeling  of  the  composite  plate  with 
patch  can  be  finished.  The  shear  test  specimens  are 
manufactured  under  the  guidelines  of  the  ASTM  D  3528  test 
procedure  as  referred  to  earlier.  The  ASTM  test  calls  for  the 
use  of  0.125  in.  and  0.063  in.  thick  adherends  for  the 
specimens.  For  this  thesis  the  coupons  were  made  of  0.190  in, 
and  0.090  in.  thick  adherends.  This  was  done  because  the 
shear  strength  of  the  adhesive  was  assumed  to  be  about  6000 
psi  and  a  test  under  these  parameters  would  have  led  to 
yielding  of  the  adherends  before  failure  of  the  adhesive. 

This  would  have  blurred  the  data  concerning  the  linear 
behavior  of  the  adhesive. 

The  surface  preparation  was  designed  with  a  cohesive 
failure  mode  in  mind.  The  adherends  are  2024-T3  bare 
aluminum.  The  surface  preparation  consists  of  a  phosphoric 
acid  anodize  followed  by  priming  with  BR-127  adhesive  primer 
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of  0.0002  in.  to  0.0003  in.  dry  film  thickness.  The  criteria 
for  a  successful  surface  preparation  was  100  percent  cohesive 
failure  in  the  coupons. 

The  adhesive  used  to  bond  the  adherends  is  3M's  AF  377 
film  adhesive.  The  film  weight  is  0.03  Ib/ft^.  The  coupons 
are  all  cut  from  the  same  bonded  plates  to  insure  uniform  cure 
conditions.  The  plates  are  cured  under  full  vacuum  (minimum 
27  in.  Hg)  in  an  oven.  The  heat-up  rate  must  be  greater  than 
40F  per  minute  up  to  205«F.  The  plates  are  held  at  that 
temperature  for  120  minutes  and  then  cooled  at  10«>F  per 
minute.  The  temperature  must  be  below  150*F  before  removing 
the  vacuum. 

The  coupons  are  cut  from  the  bonded  plates  using  a  sharp 
band  saw.  The  edges  of  the  panel  are  discarded.  After  the 
rough  cuts  are  made  the  coupons  are  milled  to  the  final 
dimensions.  A  total  of  seven  coupons  are  cut  from  the  bonded 
plates.  The  particulars  of  the  coupons  are  noted  in  Table  3. 


Specimen 

No. 

Width 

(in) 

LAP  A 
(in) 

LAP  B 
( in) 

GLT  A 
(mils) 

GLT  B 
(mils) 

4-1 

0 . 9973 

0.5045 

0.4945 

3.1 

3.0 

4-2 

0.9986 

0.5046 

0.4971 

3.0 

2.9 

4-3 

1.0023 

0.5057 

0.4990 

2 . 9 

2.8 

4-4 

0 . 9986 

0.5015 

0.4980 

3.5 

3.1 

4-5 

1.0000 

0.5018 

0.5007 

3.0 

3.1 

4-6 

0 . 9987 

0.5010 

0.5035 

3.0 

2.9 

4-7 

0 . 9960 

0.5005 

0.5060 

3.0 

3.0 

Table  3.  Double  Lap  Shear  Specimen  Dimensions 
GLT  -  Glue  Line  Thickness 
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Once  the  coupons  are  completed,  the  actual  testing 
begins.  The  test  stand  is  an  Instron  4505  (Figures  23  and 
24).  An  Instron  Strain  Gage  Extensometer  A515-11,  with  an 
initial  gage  length  of  one  inch,  is  used  to  acquire  the  strain 
data  about  the  lap  joint.  The  extensometer  spanned  the  entire 
joint.  One  knife  edge  rests  on  the  thick  adherend  and  the 
other  edge  rests  on  one  of  the  thin  adherends.  The 
extensometer  is  held  in  place  by  elastic  cord  pulling  the 
knife  edges  onto  the  specimen.  The  strain  voltage,  from  the 
extensometer,  runs  through  the  Instron  4505  before  being 
output  to  an  HP  7090A  X-Y  Plotter.  Other  output  from  the 
Instron  include,  time  of  data  acquisition,  load,  and  grip 
extension.  Only  the  load  and  strain  data  is  used  for 
evaluating  the  adhesive  properties.  The  data  is  recorded  on 
the  X-y  plotter  and  also  stored  in  an  HP  9000  Series  300 
computer  which  is  connected  to  the  plotter.  After  the  test  is 
completed  the  data  stored  in  the  computer  is  saved  onto  a 
floppy  disk. 

The  first  test  run,  4-1,  did  not  have  the  extensometer 
spanning  the  double  lap  shear  joint.  This  test  was  run  in 
order  to  find  the  approximate  failure  strength  of  the  specimen 
and  to  verify  the  material  properties  of  the  aluminum 
comprising  the  coupons.  The  extensometer  was  attached  to  the 
thick  adherend  during  this  test.  The  load  at  failure  is 
important  because  the  extensometer  can  be  damaged  if  it  is 
engaged  at  the  time  of  failure.  During  the  remaining  tests. 
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HP  7090A 
X-y  PLOTTER 


INSTRON 
STRAIN  GAGE 
EXTENSOKETER 
A515-11 
#2630-025 


LOAD 

(volts) 


STRAIN 

(volts) 


INSTRON  4505 

CROSSHEAD  SPEED  0,005  IN/MIN 


HP  9000  SERIES  300 
COMPUTER 

DATA  FRCM  X-Y  PLOTTER 
RECORDED  ONTO  FLOPPY 
DISK 


Figure  23.  Double  Lap  Shear  Test  Schematic 


Figure  24.  Double  Lap  Shear  Specimen  in  Grips 

the  extensometer  is  removed  before  the  specimen  fails.  As 
confidence  is  built  through  the  tests,  the  extensometer  is 
left  on  the  specimen  closer  and  closer  to  the  failure  load. 
The  data  acquired  reflects  this  in  that  the  data  is  complete 
for  a  percentage  of  the  load  spectrum.  As  an  example,  the 
first  test  recorded  data  until  the  load  reached  84  percent  o 
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the  failure  load,  while  the  last  test  recorded  data  until  the 
load  reached  97  percent  of  the  failure  load.  The  data  for  the 
double  lap  shear  tests  is  presented  in  Appendix  A. 

The  data  presented  graphically  is  the  lap  stress  versus 
lap  strain;  except  for  test  4-1,  which  is  normal  stress  versus 
strain  of  the  thick  adherend.  Lap  strain  is  defined  as  the 
strain  measured  along  the  length  of  the  lap  joint.  The  gaged 
portion  of  the  specimen  is  one  inch  long;  this  is  L.  The 
relative  displacement  of  the  ends  of  the  lap  joint  is 
calculated  below: 

AL  =  eL  (20) 

This  relative  displacement  is  for  two  points  along  the  edge  of 
the  specimen.  Once  the  relative  displacement,  for  an  average 
test  under  a  given  load,  is  determined,  the  CELAS2  elements 
can  be  "calibrated."  The  data  used  is  the  average  of  the 
linear  portions  of  the  stress-strain  curves  from  the  double 
lap  shear  tests.  The  spring  element  values  are  changed  until 
the  relative  displacement  of  the  two  nodes,  corresponding  to 
the  points  that  the  extensometer  knife  edges  rest  on,  matches 
the  calculated  extension. 

Panel  and  Patch  Test.  Once  the  patches  have  been  cured  and 
bonded  to  the  parent  panel,  as  described  earlier,  the  next 
step  is  the  bonding  of  strain  gages  onto  the  specimen.  The 
strain  gages  used  are  the  Micro-Measurements  Division  CEA-06- 
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125UR-350.  These  gages  are  rectangular  rosettes.  The  method 
of  attaching  the  gages  is  as  follows: 

1.  The  areas  where  gages  are  to  be  attached  are  cleaned 
thoroughly  with  acetone. 

2.  The  areas  are  then  abraded  to  the  point  where  the 
surface  is  smooth,  i.e.  free  of  fabric  imprint. 

3.  Acetone  is  again  used  to  clean  the  areas. 

4.  The  strain  gages  are  positioned  with  mylar  tape  to 
prevent  inadvertent  stretching  of  the  gages. 

5.  M-Bond  AE-10  Epoxy  is  mixed  and  applied  to  the 
underside  of  the  gages  and  to  the  area  where  the  gages 
will  be  attached. 

6.  The  gages  are  taped  back  into  place  and  pressure,  in 
the  form  of  weights,  is  applied  to  remove  excess  epoxy 
and  provide  clamp-up  pressure. 

7.  After  48  hours  the  weights  and  tape  are  removed  from 
the  gages  and  the  gages  are  cleaned  with  alcohol. 

8.  The  strain  gage  leads  are  then  soldered  to  the  tabs 
of  the  gages. 

For  the  panel  and  patch  specimens,  four  strain  gage  rosettes 
are  mounted  on  each  specimen  (for  a  total  of  12  channels  of 
strain).  One  in  the  center  of  the  patch,  two  2.25  inches  from 
the  edge  of  the  cutout,  and  one  far  field  gage,  12.5  inches 
away  from  the  cutout  edge  (Figure  25).  On  the  panels  without 
a  patch,  there  are  three  gages  mounted.  The  rosettes  are  in 
the  same  locations  except  for  the  center  one.  On  the  one 
undamaged  panel  tested,  two  gages  were  mounted  equally  spaced 
longitudinally.  All  gauges  were  attached  to  one  side  of  the 
panels,  i.e.  there  was  no  averaging  of  strain  from  the  two 
sides  of  the  specimens. 
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The  actual  dimensions  of  the  panels  tested  is  noted  in 
Table  4.  The  column  denoted  as  "DIM  A"  refers  to  the  width  of 
the  portion  of  the  panel  where  strain  gage  A  (refer  to  Figure 
25)  was  attached.  "DIM  B"  refers  to  the  width  of  the  other 
portion  of  the  panel  across  the  cutout.  Both  of  these 
dimensions  where  measured  across  the  diameter  of  the  cutout. 

If  the  panels  had  been  cut  perfectly  the  two  dimensions  would 
have  been  2.500  inches. 


Panel  No. 

DIM  B 

HOLE  DIA 

(in) 

(in) 

1-N 

N/A 

N/A 

N/A 

1-0 

2.515 

2.396 

5.003 

2-0 

2.514 

2.396 

4.999 

1-4 

2.538 

2.389 

4.995 

2-4 

2.522 

2.394 

4.998 

1-6 

2.510 

2.527 

4.994 

2-6 

2 . 512 

2.541 

4.870 

1  1-8 

2 . 527 

2.466 

5.003 

1  2-8 

2.510 

2.527 

4.994 

Table  4.  Test  Panel  Dimensions 


The  test  fixture  is  an  MTS  Testing  System  976.01-58  with 
a  Load  Cell  Model  661.31A-05  (Figures  26  and  27).  The 
specimen  is  bolted  into  adapters  for  the  hydraulic  grips.  The 
hydraulic  grips  are  about  three  inches  wide.  The  adapters 
allow  the  10  inch  wide  panels  to  be  secured  in  the  load  cell. 
The  adapter  was  originally  made  to  accept  specimens  up  to  16 
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inches  wide.  The  bolt  pattern  is  staggered  and  requires  nine 
0.5  inch  diameter  bolts  on  each  end  to  secure  the  specimen. 

The  minimum  thickness  the  adapter  will  allow  for  proper 
installation  of  a  specimen  is  0.75  inches.  Since  the  panels 
are  only  0.0832  inches  thick,  spacers  are  required  to  build  up 
the  apparent  specimen  thickness.  Half-inch  aluminum  plates 
are  used  for  this  purpose.  The  plates  are  sandblasted  and 
have  a  non-skid  surface  applied  to  the  side  in  contact  with 
the  specimen.  The  idea  behind  the  nonskid  surface  is  to  shear 
a  portion  of  the  grip  load  into  the  specimen  rather  than 
through  the  bolts.  This  way  the  bearing  loads  on  the  specimen 
can  be  reduced  without  using  time  consuming  and  costly  tabs 
which  must  be  bonded  on. 

The  strain  gage  leads  ate  i un  into  an  HP  3852A  Data 
Acquisi t 1 on/Control  Unit  {a  multiplexer).  From  the 
mu  1 t 1 p 1 fxe i  ,  the  data  is  acquired  from  an  HP  9000  Series  300 
computer  running  the  data  acquisition  software  "MULT-STR." 

Thif  data  is  displayed  in  real  time  on  the  computer  monitor. 

The  displayed  data  is  load  versus  the  various  channels  of 
stiaiii.  After  the  test  is  completed  the  data  is  saved  onto  a 
floppy  disk  foi  fuithei  evaluation.  The  data  acquisition 
.s(J^twal^•  1  .s  independent  of  the  test  cell  contiollei  and  both 
aie  started  simultaneously  at  the  beginning  of  a  test.  The 
f e;  f  ceil  contiols  the  rate  of  grip  extension  and  piovides  a 
load  ;,ignal  to  the  multiplexei  which  in  tuin  is  i  n  t  1  t  oga  t  ed 
by  the  data  acquisition  softwaie. 


STRAIN  FROM 
M-M  STRAIN 
GAGE  ROSETTE 
(volts) 


HP3S52A 

DATA  ACQUISITION/ 
CONTROL  UNIT 
(MULTIPLEXER) 


LOAD 

(volts) 


MTS  TESTING  SYSTEM 
976.01-58 

LOAD  CELL  MODEL  661.31A-05 


HP  9000  SERIES  300  COMPUTER 
DATA  ACQUISITION  SOPTWPJIE 
•MULT  STR" 

OUTPUT  RECORDED  ONTO  FLOPPY 
DISK 


Figure  26.  Panel  and  Patch  Specimen  Test  Schematic 
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The  load  cell  is  set  to  provide  grip  extension  at  0.005 
inches  per  minute.  The  scale  of  load  voltage  varies  depending 
on  the  expected  maximum  load  and  is  set  before  the  test 
begins.  The  data  acquisition  software  requires  the  following 
input  before  the  test  begins: 

1.  Test  name. 

2.  Number  of  strain  channels  to  be  received. 

3.  Gage  factor  of  the  strain  gages. 

4.  Level  of  load  per  10  volts  incoming  signal. 

5.  Scale  of  the  graphical  display  of  load  versus  strain. 

After  the  strain  data  that  has  been  collected  and  saved 

onto  a  floppy  disk,  it  is  evaluated  and  reduced  into  a  form 
comparable  to  the  data  output  from  the  FEM.  This  involves 
resolving  the  strain  gage  rosette  readings  into  Von  Wises 
strain.  Assuming  the  longitudinal  strain  gage  in  a  rosette  is 
number  1,  the  gage  45  degrees  off  from  the  longitudinal 
direction  is  number  2,  and  the  transverse  gage  is  number  3, 
the  equations  involving  these  recorded  directions  of  strain 
and  the  principal  strains  are: 

g  gmax-^-gmin^  gmax~gmin  ^^320  (21) 

^2  2 

e^=  ^  V"  COS2  )  (22) 
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(23) 


g  -  ^B>ax'*'^min 


emax-e 


“^"  cOS2  (0  +  -^) 


Strains  £2^,  62/  and  6^  are  measured  by  the  strain  gage 
rosettes.  The  principal  strains  are  calculated  after  solving 
the  above  equations  for  e^^x'  ^ndn'  Mises 

strain  is  then  obtained  from  the  principal  strains.  It's  the 
Von  Mises  strain  from  the  experiments  and  FEM  that  are 
compared  in  the  results  section  of  this  thesis.  Von  Mises 
strain  was  used  for  the  strain  field  comparison  since  this 
comparison  relates  distortional  energy  which  presents  a  more 
complete  "picture"  of  the  strain  field  behavior  at  a  given 
point.  All  tests  for  this  thesis  were  accomplished  in  the 
test  facilities  of  the  Wright  Research  and  Development  Center. 

®niax,min”  2  (24) 

«  1  1  .  -2C5+8,  , 

0=-^tan‘M— - - — -)  (25) 

2  *l”®3 


8  = 


v-^ 


2  ^  ®min 


(26) 
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VI .  Resul ts 


Finite  Element  Model.  The  results  from  the  finite  element 
model  analysis  make  intuitive  sense.  Figures  28,  29,  30,  and 
31  illustrate  the  changing  Von  Mises  strain  field  of  the 
composite  panel.  The  figures  illustrate  the  strain  field  of 
the  panel  with  the  cutout  under  load  and  the  same  panel  with 
different  patches  bonded  to  it.  The  patched  panel  figures 
show  both  surfaces  of  the  panels,  one  side  reveals  the  patch 
strain  field  entirely  while  the  other  side  shows  the  strain 
field  in  the  panel  to  the  edge  of  the  cutout. 

The  strain  fields  illustrated  reveal  how  the  load  path 
changes  from  the  unrepaired  case  to  the  repaired  case.  The  4- 
ply  patch  is  highly  stressed  due  to  the  adhesive  shearing  the 
load  into  it.  As  the  patch  laminate  begins  to  match  the 
original  panel  in  the  primary  load  direction,  the  strain 
fields  become  more  closely  matched.  The  8-ply  patched  panel 
has  only  a  slightly  higher  strain  field  than  the  far  field 
areas  of  the  panel.  Most  of  this  difference  may  be  due  to  the 
Poisson  contraction  of  the  panel  adjacent  to  the  hole  (since 
the  patches  are  all  unidirectional,  the  Poisson  contraction  is 
not  countered  jn  the  "repaired"  panels). 

The  maximum  Von  Mises  strain  in  each  of  the  four  cases 
illustrated  and  the  undamaged  case  are  noted  below: 


1 . 

Panel 

without  repair. 

^v  = 

9.04 

E-04 

2  . 

Panel 

with  4-ply 

patch. 

^v  = 

3 . 91 

E-04 

3  . 

Panel 

with  6-ply 

patch. 

^v  = 

2 . 90 

E-04 

4  . 

Panel 

with  8-ply 

patch. 

^v  = 

2.34 

E-04 

5. 

Undamaged  panel . 

^v  = 

2.05 

E-04 
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LEVELS;10  DELTA:  l.OOE-04 


Figure  29.  Strain  Field  of  Panel  with  4-Ply  Patch 
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1  .  OOE-04  1.  OOE-03 

-  ^  ^  .  --p-  )  ■)  j 

LEVELS;10  DELTA;  l.OOE-04 


Figure  30.  Strain  Field  of  Panel  with  6-Ply  Patch 
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The  areas  adjacent  to  the  cutout  reveal  the  load  tapering 
effect  of  the  adhesive.  The  strain  levels  in  the  patch  are 
raised  incrementally  from  the  edge  of  the  patch  toward  the 
edge  of  the  cutout.  This  agrees  with  intuitive  reasoning. 

Test  Specimens.  The  physical  tests  were  not  as  straight 
forward  as  the  FEM.  The  first  indication  of  problems  was  with 
the  test  of  the  "damaged"  panel  without  repair.  As  the  panel 
load  was  increased,  but  well  below  the  tensile  failure  load,  a 
portion  of  the  panel  near  the  cutout  began  to  displace  out-of- 
plane.  The  displacement  was  symmetric  about  the  longitudinal 
axis.  This  behavior  appears  to  be  similar  to  a  form  of 
transverse  buckling.  C-Scans  of  one  of  the  panels,  before 
testing,  revealed  some  delamination  around  the  cutout.  The 
delaminated  areas  extended  into  the  panel  no  more  than  1/8 
inch.  The  delamination  was  probably  caused  in  the  machining 
of  the  cutout.  Under  the  transverse  compressive  Poisson 
stress  near  the  cutout,  the  delamination  may  have  caused 
transverse  shear  which  may  have  led  to  the  transverse 
instability.  This  instability  may  also  have  been  aggravated 
by  the  relatively  thin  panel  and  large  cutout  compared  with 
the  other  panel  dimensions.  This  is  indicated  in  Figures  34 
and  35.  With  the  load  between  8000  and  9000  pounds  the  strain 
at  point  B  changed  sharply  (refer  to  Figure  25).  This 
behavior  was  not  modeled  and  it  limited  the  acceptable  data, 
of  the  test,  to  the  moments  before  the  out-of-plane  behavior 
became  the  dominant  strain  mechanism.  Generally,  this  limited 
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the  useful  data  to  end  loads  below  10,000  pounds.  While  the 
unrepaired  panel  was  the  most  affected  by  this  phenomenon,  all 
the  panels  revealed  signs  of  "transverse  buckling"  during  the 
tests.  The  typical  failure  mode  of  the  patched  panels  was  one 
of  transverse  bending  failure  of  the  patches  accompanied  by 
delamination  of  the  patch  from  the  panel  and  then  tensile 
failure,  aggravated  by  out-of-plane  displacement,  of  the 
parent  panel.  Since  the  patch  fibers  were  all  aligned  with 
the  load  axis,  the  patch  transverse  strength  was  basically 
limited  to  that  of  the  matrix.  Figure  32  illustrates  the 
patch  transverse  failure  mode. 

During  the  testing,  all  the  patches  underwent  local 
bending  near  the  edge  of  the  cutout.  The  patches  "sank"  into 
the  cutout  after  the  load  was  applied  to  the  panel.  This 
bending  did  not  appear  to  affect  the  centrally  located  strain 
gage  rosette. 

All  the  strain  gage  rosettes  remained  attached  during  the 
tests  and  acceptable  data  was  generated  through  them.  The  Von 
Mises  strain  at  the  gages  was  measured  and  correlated  with  the 
predicted  strain  (Figures  33-41).  The  gages  attached  at 
points  A  and  B  (refer  to  Figure  25)  on  the  patched  panels  did 
not  correlate  well  with  the  predicted  strain.  This  is 
probably  due  to  the  complex  three-dimensional  load  interaction 
of  the  patch  and  panel  at  those  locations,  i.e.  the  local 
patch  bending  and  load  transfer.  The  data  presented  for  the 
patched  panels  reflects  that  taken  from  gages  centered  on  the 
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patch  and  far  field,  points  D  and  C  (refer  to  Figure  25). 

Figures  36-41  show  clearly  the  excellent  correlation 
between  the  analytical  model  and  experimental  data.  The  use 
of  the  CELAS2  elements  in  this  elastically  loaded  situation 
clearly  modeled  the  elastic  behavior  of  the  panel -adhesive- 
patch  load  transfer  mechanism  successfully.  The  slight 
variation  among  the  curves  may  be  due  to  the  imperfect  test 
specimens  as  noted  in  Table  4. 
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Figure  33.  Von  Mises  Strain  versus  Load,  Undamaged  Panel 
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Figure  34.  Von  Mises  Strain  versus  Load,  Unrepaired  Panel 


Number  One 
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Analytic 


Experimental 


L 


35.  Von  Mises  Strain  versus  Load,  Unrepaired  Panel 

Number  Two 
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STRA  : 


Figure  37.  Von  Mises  Strain  versus  Load,  4-Ply  Patch  Panel 

Number  Two 
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Figure  38.  Von  Mises  Strain  versus  Load,  6-Ply  Patch  Panel 

Number  One 
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igure  39.  Von  Mises  Strain  versus  Load,  6-Ply  Patch  Panel 


Number  Two 
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Analytic 


Experiment 


Figure  40.  Von  Mises  Strain  versus  Load,  8-Ply  Patch  Panel 

Number  One 
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Figure  41.  Von  Mises  Strain  versus  Load,  8-Ply  Patch  Panel 

Number  Two 
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VII.  Conclusions 


The  FEM .  This  thesis  has  shown  conclusively  that  using  CELAS2 
spring  elements  to  model  an  adhesively  bonded  joint  can 
adequately  predict  a  basically  two-dimensional  load  transfer 
process.  Both  orthotropic  and  isotropic  materials  were  used 
and  modeled  in  this  thesis.  Use  of  the  MSC/NASTRAN  CELAS2 
spring  elements  proved  satisfactory  in  modeling  the  adhesive 
bond  in  both  cases.  The  FEM  accurately  predicted  the  in¬ 
plane  behavior  of  the  damaged  and  repaired  specimens.  Good 
correlation  was  observed  both  in  displacement  and  strain. 

The  effect  of  the  bond  was  accurately  modeled  through  the 
use  of  the  linear  spring  elements.  The  use  of  linear  spring 
elements,  in  an  FEM,  provides  insight  into  the  load  transfer 
occurring  in  a  bonded  joint.  This  visualization  should  be 
extremely  helpful  is  analyzing  complex  geometry  bonded  joints. 
The  modeling  method  itself  is  not  exceptionally  difficult  or 
time  consuming  and  yet  it  provides  realistic  answers  to 
complex  problems. 

This  thesis  did  not  explore  the  three-dimensional 
possibilities  of  the  CELAS2  element.  Since  this  element  has 
three-dimensional  capabilities,  it  may  be  feasible  and 
realistic,  to  model  such  cases  as  a  thin  adherend  single  lap 
shear  test  specimen  which  undergoes  shear  and  peel  during  the 
test.  These  specimens  generally  fail  through  the  combined 
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action  of  in-plane  and  out-of-plane  stress.  3uch  a  modeling 
and  test  analysis  may  be  fairly  straight  forward. 

MSC/NASTRAN  is  also  equipped  with  rotational  spring 
elements  for  modeling  purposes.  The  spring  force  generated  in 
the  element  is  dependent  on  the  rotation  of  the  node  it  is 
associated  with.  The  element  does  not  function  between 
elements  but  between  particular  nodes  and  the  global 
coordinate  system.  The  rotational  spring  applies  a  moment  to 
the  node.  If  the  rotational  spring  element  depended  on  the 
relative  rotation  of  two  associated  nodes,  then  the  shear 
modulus  of  the  adhesive  could  have  been  used  directly  to  model 
the  bond. 

The  Materials.  To  be  fair,  nothing  good  or  bad  can  be 
concluded  about  the  low  energy  prepreg  SP  377  CUl  material 
provided.  At  the  time  of  this  writing,  3M  is  still  not  sure 
how  the  defective  roll  of  prepreg  managed  to  leave  their 
production  facility.  Other  organizations,  namely  the  U.S. 
Navy,  have  received  and  tested  this  material.  The  Navy 
reported  some  differences  between  the  published  SP  377  CUl 
mechanical  properties  and  what  they  observed  in  their  tests. 
Further  research  into  the  applicability  of  this  material  for 
ABDR  is  necessary  before  any  conclusions  one  way  or  the  other 
may  be  made. 

The  AF  377  adhesive  film  worked  well  during  the  double 
lap  shear  specimen  testing  and  the  composite  panel  and  patch 
fabrication  and  testing.  The  adhesive  provided  an  easily 
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controlled  and  curable  bonding  agent.  The  adhesive  shear 
strength  data  in  Appendix  A  shows  the  impressive  capabilities 
of  this  resin  system. 

The  hot  bonder  concept  was  tried  but  no  conclusive 
arguments  can  be  made  either  way  due  to  the  defective  raw 
materials  involved. 
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Appendix  A:  Double  Lap  Shear  Test  Results  Data 


The  plate  thickness  of  the  thick  and  thin  adherends  is 
0.1965  in.  and  0.0902  in.  respectively. 


Specimen 
No . 

' 

Failure 

Load 

(lbs) 

Thick 

Adherend 

Stress 

(psi) 

Thin 

Adherend 

Stress 

(psi) 

Adhesive 

Shear 

L'tress 

(psi) 

4-1 

6126 

31260 

34050 

6149 

6086 

31020 

33780 

6084 

mSm 

6110 

31020 

33790 

6067 

4-4 

6068 

30920 

33680 

6080 

4-5 

6150 

31300 

34090 

6135 

4-6 

31060 

33830 

6076 

4-7 

6058 

30950 

33720 

6043 

Table  5.  Double  Lap  Shear  Failure  Data 


Figure  42.  Test  4-1  Thick  Adherend  Normal  Stress  Versus 

Normal  Strain 
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Figure  48.  Test  4-7  Shear  Stress  Versus  Lap  Strain 
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